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Abstract 

Using the recently proposed ggr]' effective vertex, we investigate the pro- 
duction of rj' from gluon fusion in pA colhsions. We show that measuring 
r]' production cross-section at moderate x^' yields direct information on the 
small X gluon distribution function of the nucleus. At RHIC, the smallest 
accessible x turns out to be O(10~^) and at LHC, it is 0(10^^). Therefore, 
rj' is an excellent probe of the Color Glass Condensate. 



1 Introduction 



The meson has many interesting properties, the most discussed being its unusual 
mass. By spontaneously breaking ?7(3) symmetry of the hadronic Lagrangian, one 
would expect that r( would be a Goldstone boson whose finite mass is due to the 
finite mass of the quarks. In reality, the rj mass (M = 0.958 GeV) is much heavier 
than the mass of the t] {M^j = 0.547 GeV). The quark model Gell-Mann-Okubo 
formula for the singlet and the octet mass 

Mo = T^H-f+^o+m^^) (1) 

and 

M| = - {2ml+ + 2mlo + 2ml+ + mlo) (2) 

is incapable of explaining this large mass difference between the (mostly) flavor octet 
state 7] and the (mostly) singlet state rj'. In fact, the singlet mass Mq is smaller than 
the octet mass Mg. Therefore, in the zero quark mass limit, about half of the r]' 
mass remains while the octet pseudo-scalar meson (vr, K, rj) masses go to zero. 

The resolution of this Ua{X) problem was provided by 't Hooft and Witten 
p|. In these works, it was argued that the anomalously large mass of rj' is due to 
the mixing of the gluon state with the flavor-singlet quark state through the axial 
anomaly 

d^Jl = 2iY^ rufqp.qf + 2iV;-|^ Tr {G,^G^^) (3) 
/ 

with non- vanishing (GG^. This implies that the r]' mesons have a large glue content. 
It also implies that the gluon fusion process gg —>■ rj' is possible. 

In this paper, we study the generation of rj' mesons from such a gluon fusion 
process in proton- nucleus '(pA)' collisions. In particular, we argue that measuring 
the rj' meson momentum spectrum enables us to have a direct access to the gluon 
density in the small x region of the heavy nucleus. This is made possible by using a 
recently proposed effective vertex by Atwood and Soni |^ between gluons and the 
7]' meson. 

In 1^, as a step to explain an anomalously large B to rj' branching fraction, the 
authors wrote down a Wess-Zumino- Witten type of effective vertex: 

T,,5"^ = H{p', q', P') 6^' (e^) a (ej), (4) 
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roughly corresponding to an interaction Lagrangian of the form L ~ rj'GG. Here, 
p, q are the gluons momenta, e^'^ are the associated polarization vectors, and the 
jafe (jgnotes the singlet combination of the gluons. From now on, we denote the rj' 
momentum by a capital letter P. The factor HijP', g^, P^) is a form factor. By 
analyzing J/ip ^ ri''~^ process, the authors of [§] then estimate that in the on-shell 
limit 

Hq = H{0, 0, A/P) ^ 1.8 GeV"^ (5) 

where M = 0.958 GeV is the mass of the rj' meson. This vertex was then used to 
analyze the 5-boson decay to 1]' assuming that the and dependence of H is 
weakQ. This form of effective vertex has been also confirmed by a calculation of the 
triangle diagram contribution to the axial anomaly 0. 

The above effective vertex represents a hadronization matrix element between 
a two-gluon state and a hadron state. To the authors' knowledge, this is unique. 
Certainly, many constituent quark models can relate valence quarks to hadrons. 
However, except for the above ggrj' effective vertex, there is no other known matrix 
element between the gluons and a known hadron state. We also note that the 
momenta involved in this process is not so soft compared to Aqcd or the pion mass 
since the rj' mass is almost 1 GeV. 

Recently, these interesting features of the ggr]' vertex were exploited to calculate 
the 7]' production from a hadronizing quark-gluon plasma 0. In that study, a 
Boltzmann equation was used to evolve the rj' density towards the hadronization 
time in heavy ion collisions. To use the ggr]' vertex in this 'AA' environment, a 
certain set of assumptions had to be made. These assumptions mainly involve the 
as-yet-unknown in-medium properties of the particles and vertices involved in the 
calculation. 

In this paper, we focus on production of rj' in pA collisions for two reasons. First, 
much is known about the structure of proton from DIS experiments at HERA and 
elsewhere [0. Assuming that we understand the parton distributions in a proton, 
pA experiments are an excellent tool with which to investigate parton distributions 
of nuclei at small x. This becomes even more important in light of the fact that 
not much is known about modification of parton distributions in nuclei at small 

^ This weak dependence assumption has been questioned by Kagan and Petrov who argued 
that when one of the gluon momenta is ofF-shell, the form factor has to be replaced by 

Hiq',p', M') « HoMyiM^ - q') . (6) 

However, this is not relevant to the present study. 
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X and intermediate since all existing fixed target eA experiments have limited 
coverage in x and p, 0, |Tl|, |12|. This region of small x and intermediate 



1 — 10 GeV^) is the domain of high gluon density QCD and semi-hard physics 
where interesting phenomena such as saturation of gluons, etc. are expected to 
occur. pA experiments will give us a new venue, in addition to eA and AA, in which 
to investigate the high gluon density phase of QCD. 

Also, pA collisions are much cleaner than A A in the sense that one can avoid 
many of the complications due to the in-medium effects. We do not expect to have 
a Quark Gluon Plasma (or a hot hadronic phase) created in pA collisions and can 
therefore avoid the difficulties associated with understanding the possible formation 
of plasma and its detailed properties. In this work, we will show how one can extract 
the nuclear gluon distribution function at small x from pA experiments at RHIC 
assuming one knows the gluon distribution function in a proton at not too small 
values (0(a; ~ 0.1 - 0.01)). 

It should be noted that even though our experimental knowledge of the parton 
distribution functions in a proton is good, due mostly to HERA, in some pp induced 
processes such as pp —>■ jX or pp — > nX, it is quite common to introduce an 
"intrinsic" momentum (k'f) into the standard parton distribution functions in order 



to fit the experimental data |T3[ . For example, 

xG{x,Q^) = J d^ktxG{x,Q\k^) = J d^ktxG{x,Q'^)f{k'^) (7) 
where f{k^) is usually taken to be a Gaussian with width (k^). More explicitly, 

so that / d'^ktf{k^) = 1. The value of (fc^) is extracted from best fits to experimental 
data on dileptons, dijets, etc. and can be as large as 1 — 5 GeV^ at the Tevatron. 
Its value depends also on the energy of the collision and the process considered and 
is expected to be smaller for production of heavier particles. 

The origin of this "intrinsic momentum" is mostly initial state radiation of soft 
gluons. A rigorous investigation of these soft emissions has been carried out in some 
processes, such as dilepton production, which effectively introduces this "intrinsic" 
momentum into the standard coUinear factorization based cross sections WM. In 



addition to the "intrinsic" momentum due to the initial state radiation, there is a 
genuine non-perturbative contribution which is energy and process independent and 
is of order 0(200 — 300 MeV) consistent with the uncertainty principle. Phenomeno- 
logically, these results can also be obtained by using the modified parton distribution 
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functions as defined in (|^) without the resummation of the soft gluon radiation. This 
is especially useful since this soft gluon resummation has been carried out only for 
a limited number of processes. Nevertheless, this resummation is usually taken to 
be the theoretical justification for introduction of this "intrinsic" momentum into 
the standard parton distribution functions. There is also the Cronin effect in pA 
collisions which can contribute to this intrinsic momentum. Even though Cronin 
effect is not well understood in pQCD, one can phenomenologically understand it 
in terms of multiple scattering of the proton on the nucleus which causes the pt 



broadening of the proton [|15 



In large nuclei and at small x, one also has to include high gluon density effects. 



In |16, 17, p8l O, EO] it is shown that high gluon densities lead to a potentially large 



intrinsic scale Ql{x) ~ A^^^ ^.2-.4 in large nuclei. Most gluons in the wave function 
of a nucleus in a high energy collision have momenta of the order of this scale. In 
this paper we will investigate production of 7]' in pA collisions as means of extracting 
the gluon distribution function of nuclei at small x and study its dependence on the 
intrinsic momentum (kf) in protons and on saturation scale Ql{x) in nuclei. 



2 r]' Production Cross Section 

In this section, we will calculate the t]' production in pA collisions and show how 
to extract the gluon distribution function of the nucleus A. We will consider three 
distinct cases; first, we will calculate the rj' production cross section, with the stan- 
dard gluon distribution function, and using the coUinear factorization formalism. 
Then, we will consider the case when the gluons in the nucleus A are at small x 
and include saturation effects in the nuclear gluon distribution function. Finally, 
we will allow for "intrinsic" momentum in a proton and consider the case when the 
gluon distribution function of the proton is also modified to take the "intrinsic" 
momentum of the proton into account. 

2.1 Collinear Factorized Cross Section 

Using collinear factorization theorems in perturbative QCD, one can write the t]' 
production cross section in pA collisions as 

daP^^v'x ^ J Qp^^^^ q2^ G'^{x2, Q)) da^^^"^' (9) 

where G'(x2,Qj) are the gluon distribution function of the proton and 

nucleus respectively and Xi and X2 are the momentum fraction of the incoming 
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gluons while Xl = 2Pz/^/s and Xe = 2Epj^ are the momentum and energy 
fractions of the produced r( respectively. The distribution functions, G(xi, Q^) and 
G(x2, Q/); depend on a factorization scale which is shown explicitly and will be 
taken to be ~ M^, in numerical calculations. We also note that there are no intrinsic 
momenta included in either distribution functions since it is not, strictly speaking, 
allowed in collinear factorization formalism. 

The partonic cross-section for gg rj' is given by 

= 777^7^;^|i^l^^^-'''l'^2^^^^^ (10) 

where p and q are the momenta of the incoming gluons while P is the momentum of 
the produced rj'. The matrix element squared, after averaging over the initial spin 
and color degrees of freedom, becomes 

\T„'f = j^H',M'^, (11) 



Using (11), the differential cross section becomes 



This is the elementary partonic cross section which goes into our calculation of 
7]' production in pA collisions. It should be noted that collinear factorization for 
7]' production has not been explicitly proven but is analogous to using (rigorously 
proven) collinear factorization in production of heavy quarks, high pt jets, high 
mass dileptons, etc. and that since the natural scale involved [rj' mass) is not very 
large, we may have large higher order corrections. This reflects itself in having cross 
sections which are quite sensitive to a change of factorization scale. Therefore, we 
feel our results are most reliable for the ratio of cross sections {pA over pp) rather 
than the absolute cross sections. 

Using Eq. ([12|) and resolving the delta functions gives 
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° x+GP(x+,g5)x_G^(x_,g5) (13) 



where 



64 Xe 



xe±xl Ep±P^ 

x± = 7> = ^ 14 

2 
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In terms of the rj' rapidity y or the longitudinal momentum fraction xl, these are 



M„, e^y + /s ± XL 
= = 2 ^''^ 

Here we take the positive y (as well as xl) to correspond to the proton direction. 
For xl 3> Mrj' I \fs^ we can also write 

x+ ~ Xl and a;_ ~ — ^ (16) 

SXl 

Equation ( PTSD is a very simple and useful relation which enables us to extract the 
gluon distribution function of the nucleus directly (without any de-convolution, as 
is usually the case) from the experimentally measured Vj production cross section. 
Here we are assuming that the gluon distribution function of a proton is known from 
HERA to good accuracy |^ (to better than 10%) which is the case unless x+ > 0.2. 
Also, in the region x+ < 10"'', one also has a large uncertainty in determining the 
gluon distribution function in a proton but we avoid this region by considering small 
values of x_ (recall x_ is the momentum fraction of gluons in the nucleus) since this 
is the region where nuclear modification of the gluon distribution function is most 
pronounced and least known. 

Let us make a rough estimate of the x values one can probe in a high energy 
collision, such as those at RHIC or LHC Since there are no experimental data on 
nuclear gluon distribution function at x < 10^^ for > \ GeV, we would like to 
focus on the small x (in the nucleus) kinematic region. This corresponds to the case 
when X- in (|l^) is small. This in turn corresponds to the case when the measured 
Vj has a large (much larger than its mass or Pj) longitudinal momentum, i.e. in the 
forward region (c./. Eq.(p!5|)). 

At RHIC, the largest r( rapidity in the CM frame is about 5.36. Therefore in 
principle, the x_ range is 

2.2 X 10"^ <X-< 0.005 (17) 

using Eq. (p!5|) . Since = MVs, this range of x_ corresponds to 

0.0005 < x+ < 1.0 (18) 

For instance, if we detect an rj' at y = 3.04, that would correspond to x^ = 0.10 
and X- = 2.3 x 10~^. At LHC, the largest t]' rapidity in the CM frame is about 8.6. 
This implies that 

3.2 X IQ-^ <x^< 1.7 X 10-^ (19) 
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Since = MVs, this range of x_ corresponds to 

1.7 X 10"^ < x+ < 1.0 (20) 

Again, if we fix a;+ = 0.1 or y^/ = 6.35, then x_ = 3.0 x 10^^. These would be, by 
orders of magnitude, the smallest values of x where gluon distribution in a nucleus 
has ever been measured. For example, the smallest values of x measurable in fixed 
target DIS experiments at CERN and Fermilab is x ~ 10~^ at similar [22, 23 . 
Therefore, by measuring the r/' production cross section at large rapidities, one can, 
for the first time, determine the nuclear gluon distribution function at very small 
X and moderately large Q^. This will help one understand the nature of nuclear 
shadowing in QCD and determine the role of high parton density (higher twist) 
effects in nuclear shadowing. 

2.2 Intrinsic Momenta 

In this section, we will include the effects of intrinsic momenta in the gluon distri- 
bution of both proton and nucleus and investigate the dependence of t]' production 
cross section on the intrinsic momenta. Since introducing the intrinsic momentum 
is done in the phenomenology spirit as discussed earlier, we will take its width to 
be the saturation scale, Qs in nucleus. This comes about because at small values of 
X in nuclei, one encounters the high gluon density region of QCD where non-linear 
gluon recombination effects become important |25|, ^ |26| , |28| , The high 



gluon density introduces a new scale, the saturation scale Qs{x), which grows with 
energy {Q^x) > Aqcd at small x) |0|, |T], |2|, |3|, ||, ||. Most gluons 



in the wave function of a nucleus have momenta of the order of Qs and therefore, 
one may take the average (the width of the Gaussian) intrinsic momentum to be 
this saturation scale Qs- The value of Q^ is estimated to be ~ 1 — 2 GeV^ at RHIC. 
As a first approximation and to keep our expressions simple, we will ignore the x 
dependence of Q^ and calculate the r]' production cross section with different values 
of Qs- 

As a first step, we will include the intrinsic momentum in a nucleus and not in 
a proton since the nuclear intrinsic momentum is expected to be much larger than 
that of a proton. Later, we will consider the most general case. The momenta of the 
incoming gluons are now, including the intrinsic momentum pt ~ Qs in the nucleus, 

g= (xiv^/2,0t,Xiv^/2) (21) 

and 

P = (x2v^/2 + p,V2x2v^,Pt, -X2Vs/2 + p'f/2x2Vs) (22) 
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The r( production cross section then becomes 



^^vA^^x ^ J ^^^^^^ ^2^^ Qp(^^^^ G^{x2, Q},Pt) da"' (23) 
where G^{x,Qj,pf) = G'^{x,Q'^) f{pf) and fipf) is a Gaussian (|) with width of 

Ql 

Substituting Eq. (|T2|) and using the 6 functions gives 

,,p , = ^x.C'ixu Q})x,G\x2, Q}, Pt) (24) 

where xi = M^,/sx- and X2 = X- and Pt is the measured transverse momentum 
of the produced rj'. Therefore, measuring the produced t]' at different rapidities 
{xl) and transverse momentum would determine the gluon distribution of nucleus 
at various x and scale Q/ ~ M^i. 

To get a feeling for the kinematics region where one can extract the nuclear gluon 
distribution function, let us set Pt = 1 GeV for RHIC. The range of accessible X2 
for the nucleus is now 

3.3 X 10"^ <X2< 0.007 (RHIC) (25) 

which is still very small. For LHC, we have access to higher transverse moment. If 
we set Pt = 5 GeV. The range of X2 in this case is 

1.7 X 10"^ <X2< 9.3 X 10"^ (LHC) (26) 

Even for this large Pt, the accessible X2 range goes much beyond the currently 
accessible region. 

Let us now consider the most general case when both the incoming proton and 
nucleus have intrinsic momenta. It is still expected that the intrinsic momentum 
of the proton should be less than that of a nucleus. In principle, we could also 
introduce a saturation scale for protons since at high enough energies (small x), 
the non-linearities of the gluonic fields which lead to the saturation scale would 
become important even in a proton. However, in this work, we would like to restrict 
ourselves to not too small values of a; in a proton so that we can explore the very 
small X region in a nucleus. Therefore, we will ignore high gluon densities effects in 
a proton and take the average intrinsic momentum in the proton as determined in 
pp experiments. 
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In the case when both proton and nucleus have transverse momenta, the mo- 
menta of the incoming gluons become 



q = + ql/2xi^s, qt, x^^sjl - ql/2xi^s) (27) 

and 

P = (x2v^/2 + p?/2x2v^, Pu -X2^s/2 + pI/2x2^s) (28) 
The differential cross section then can be written as 

^^PAWX _ j ^^^^^^ ^2p^ ^2^^ Qp^^^^ q2^ ^A^^^^ q2^ (39) 

which, upon using Eq. (^), becomes 

dxL(PPt 64 xe 

where qt = Pt — Pt and 

1 



/ Spt x^G^xi, QJ, qt) x,G\x2, Q}, Pt) fr-T-pO) 



M\pt, qt) = -[M^ + 2pfqt + ^{M^ + 2pfqt)' - Ap^i ] (31) 



with 



Xi = Xo = — , 32 

sx_ M^ + qf ^ ' 

using x± = {xe ± xl)/2. The pt integration can be performed numerically which 
would then directly relate the gluon distribution function in a nucleus G^{x, Q'j) at 
X and ~ M^, to the measured rapidity and Pt of the produced rj'. 

However, since the t]' production cross section in pA collisions is not experimen- 
tally known, we can conversely use the available parameterizations of the nuclear 
gluon distribution function to predict the t]' production cross section in pA collisions, 
for example, in those planned at RHIC. In this work, we will take this approach and 
use two available parameterizations of the gluon distribution function in nuclei due 



to 111 and 12 



In Figure (0) we show the double differential cross section ^^^^2p^ as a function 
of 1]' transverse momentum at xl = 0.1 for two different values of the intrinsic 
momentum (saturation scale) in the nucleus. As is seen, increasing the intrinsic 
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Figure 1: Vj production cross section in "pA at Xi^ = 0.1 and y/s = 200 vs. rj' 
transverse momentum P^. 



momentum reduces the differential cross section by a factor of 3 — 4. It sliould be 
noted that were we to plot our results vs. rapidity, the cross section would increase 
with increasing intrinsic momentum. 

In Figure we show the Pt integrated cross section vs. the rj' longitudinal 
momentum ratio x^. Again a decrease in the cross section is seen as we go to higher 
intrinsic momenta. In both Figures (|I]) and BQV shadowing is used. 

In Figure we show the ratio of rj' production cross sections in pp and pA 
collisions for two representative values of the nuclear intrinsic momenta. As shown, 
this ratio is quite sensitive to the value of the intrinsic momenta. We have checked 
that the choice of parameterization of nuclear shadowing (||ll| or |jl2|) makes very 
little difference. 

Following the thermal model estimate, we expect that the total cross-section for 
Tj' production should be about 3-5% of 7r° cross-section. Among the decay mode of 
rj', the experimentally cleanest one is the rj' — >■ 77 channel. This process has the 
branching ratio of about 2 %. Therefore, the branching ratio times the cross-section 
should be on the order of Aa-j^o/lOOO ~ Imb. To measure this at RHIC, one will 
have to separate these photons from those coming from other sources. Nevertheless, 
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Figure 2: rj' production cross section at ^/s = 200 vs. xl- 





1.00 r 




0.90 - 




0.80 - 




0.70 - 


X 


■o 




■o 


0.60 - 


< 






0.50 - 






■o 






c 


0.40 - 








0.30 - 




0.20 - 




0.10 - 



BQV shadowing and p ,=1 .0 GeV 
BQV sliadowing and pV2 GeV 



0.001 




0.01 



0.1 



X, 



Figure 3: Ratio of pA to pp rj' production cross sections at i/i = 200. 
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both STAR and PHENIX detectors will be able to detect ?7's. As a very rough 
estimate, PHENIX may be able to measure rj's in the transverse momentum range 
1 < Pt < 5GeV and (pseudo)rapidity range —0.3 to 0.3p9[]. 



3 Conclusion 

We have investigated the production of r]' in pA collisions and its sensitivity to the 
intrinsic momenta (saturation scale) in nuclei. We have shown that experimental 
measurement of this cross section at RHIC would lead to determination of the 
gluon distribution function in nuclei in a much wider kinematic region accessible 
by any other process in the current experiments. Conversely, we have used the 
available parameterizations of the gluon distribution function in nuclei to predict 
the 7]' production cross section in pA collisions at RHIC. 

There are several points which need to be better understood in order to make 
our calculation more accurate. First of all, we have included only the gluon fusion 
and ignored all other possible mechanisms of i]' production. For example, it may be 
possible to produce rj' through production of quarks and anti-quarks which would 
hadronize into a r]' through a completely non-perturbative fragmentation function. 
The effect of this mechanism is hard to estimate since quark (or anti-quark) frag- 
mentation into 7]' has not been measured. 

We have also taken the r]' mass to be a constant (958 MeV) since we don't expect 
to have a Quark Gluon Plasma produced in a pp or pA collision. One then would 
have had to investigate the temperature dependence of the anomaly and include 
the poorly understood temperature dependence of the t]' mass. However, there 
is a recent study of high gluon density effects on instantons which indicates 
the average instanton size shrinks as the saturation scale Qs grows analogous to 
the finite temperature case. This would imply that t]' would become lighter and 
its production cross section would increase at higher energies due to the rise of Qs- 
Decreasing t]' mass would reduce our ggrj' matrix element (0) and our rj' production 
cross section. Therefore measuring rj' production cross section at different energies, 
for example at RHIC would help clarify this point. 

Also, we are using the collinear factorization for a process for which there is no 
proof even though one may naively expect it to hold in analogy to heavy quark pro- 
duction. Our factorization scale, M^, is not very high and our results are somewhat 
sensitive to the change in this scale. Therefore, we feel that our results are most 
reliable for ratio of pA to pp cross sections for rj' production. 
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